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Abstract 


This  report,  as  the  title  indicates,  is  an  extension  of  the  work 
begun  in  Research  Report  No.  WP-1, 

The  approximation  method  described  in  WP-1  is  applied  here  to  the 
simple  problem  of  reflection  from  a  plane  interface  separating  two  homogeneous 
media  with  different  dielectric  constants.  The  resulting  approximate  reflection 
coefficient  is  the  first  term  of  a  rapidly  convergent  series  whose  sum  is 
shown  to  be  the  standard  exact  expression  for  the  reflection  coefficient  of 
this  problem. 

In  the  second  part  of  the  present  work  a  different  approximation 
method  is  described.  The  new  method  is  appropriate  for  certain  cases  in  which 
the  previous  method  fails  or  gives  poor  results,  e.g.,  the  esse  of  a  wave 
incident  at  a  nearly  grazing  angle. 
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1.  Introduction 

In  Research  Report  WP-1  a  method  for  obtaining  the  approximate  reflection 
coefficient  of  an  inhomogeneous  transmission  line  was  described.  An  estimate 
of  the  error  produced  by  the  approximation  was  given,  and  it  was  seen  that  the 
error  depended  only  on  a  quantity  which  was  called  the  total  variation  of  the 
characteristic  impedance  throughout  the  medium.  Thus,  we  refer  to  the  approxi- 
mation of  report  WP-1  as  the  small  variation  approximation. 

A  basic  assumption  in  the  derivation  of  the  small  variation  approximation 
was  that  the  characteristic  impedance  be  a  continuous  function  of  distance  along 
the  line.  Since  the  error  estimate  apparently  has  nothing  to  do  with  the  con- 
tinuity of  the  characteristic  impedance,  we  might  expect  that  the  method  can 
be  extended  to  cover  a  more  general  class  of  problems.  Indeed,  we  shall  demon- 
strate in  Section  2  of  this  report  that  the  characteristic  impedance  is  allowed 
to  have  i\mp  discontinuities,  and  we  shall  indicate  how  the  result  of  report 
WP-1  can  be  extended  to  cover  the  case  of  an  arbitrary  medium  whose  characteristic 
impedance  is  piecewise  continuous.  A  simple  example  will  be  treated  in  detail 
to  illustrate  the  application  of  the  small  variation  approximation  as  well  as 
the  significance  of  our  error  estimate. 

In  Section  3  we  shall  describe  another  class  of  approximations  for  the 
wave  solutions  of  an  inhomogeneous  transmission  line.  The  new  class  of  approxi- 
mations will  apply  in  some  cases  to  which  the  small  variation  approximation 
cannot  be  applied:  in  particular,  the  case  of  a  plane  wave  incident  at  an  angle 
near  ninety  degrees  on  an  inhomogeneous  plane  stratified  medium.  An  example 
of  the  new  method  applied  to  this  problem  will  be  given  in  Section  h» 
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It  can  be  seen  from  this  description  that  Section  2  is  independent  of 
Sections  3  and  U»  except  insofar  as  all  three  sections  refer  to  the  results 
of  report  WP-1. 

2,     Example  of  the  small  variation  perturbation 

In  order  to  show  in  detail  the  exact  nature  of  the  small  variation 
approximation  we  shall  demonstrate  its  application  to  a  very  simple  problem 
whose  solution  is  well  known.     The  example  we  shall  use  devolves  upon  the  case 
of  two  different  homogeneous  transmission  lines   joined  at  the  point  of  dis- 
continuity z  «  0,   the  line  in  the  region  z  <  0  having  the  constant  characteristic 
impedance  Ji_   and  the  line  in  the  region  z  >  0  having  the  constant  characteristic 
impedance  S^p. 

The  well-known  exact  relation  for  the  reflection  coefficient  of  this 
pair  of  transmission  lines  at  the  point  z  =  0  is 

(2.1)  R  =  (^2-^l^/^^2*^l^* 

It  is  possible  to  derive  (2,1)  from  the  integral  equations  of  report  WP-l; 

/A  (^ 

(2.2)  u(z)  =  u(A)  +  [  {5'(y)/2^(y)]  exp^2jjk(8)d3|d(y)dy 


y 

'A 


d(z)  - 


j  (5^'(y)/k(y)|  exp|2Wk(s)ds}u(y)dy 


for  the  upgoing  and  downgoing  factors  u(z)  and  d(z)  of  the  traveling  waves j 
the  reflection  coefficient  at  the  point  z  =  0  is  given  in  terms  of  these 
functions  by 
(2.3)  R  =  d(O)MO). 
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In  report  WP-1  we  derived  equations  (2,2)  under  the  assumption  that 
the  characteristic  impedance  J^(z)  is  a  continuous  function,  whereas  in  the 
present  case  ^(z)  has  a  jump  discontinuity  at  z  «=  0.  Thus,  it  will  be  necessary 
to  treat  the  transmission  line  of  our  example  as  the  limiting  form  of  a  sequence 
of  continuous  transmission  lines.  This  procedure  is  also  obviously  necessitated 
by  the  presence  of  the  derivative  of  ^(z)  in  both  integrals  on  the  right  of 
(2,2).  Each  one  of  the  continuous  transmission  lines  in  our  assumed  sequence 
is  homogeneous  except  in  an  interval  0  <  z  <  6  where  }C(z)  varies  continuously 
from  J^-  to  t^^.     In  the  limit  6-^0  the  sequence  converges  to  the  case  of  our 
example. 

If  we  carry  out  the  infinite  iteration  of  (2.2)  described  in  report  WP-1, 
the  iteration  whose  first  order  term  constitutes  the  approximate  solution  of  the 
small  variation  method,  we  shall  obtain  the  functions  u(z)  and  d(z)  exactly 
and  thus  from  (2»3)  the  exact  reflection  coefficient  at  z  =0.  The  infinite 
iteration  is  too  complicated  to  apply  exactly  to  any  one  of  the  sequences  of 
continuous  transmission  line  problems.  However,  we  can  obtain  an  estimate 
fiT>m  which  the  limiting  form  of  the  sequence  of  solutions  corresponding  to  the 
sequence  of  continuous  transmission  lines  can  be  obtained.  The  limit  of  the 
resulting  sequence  of  reflection  coefficients  will  turn  out  to  be  identical 
with  the  reflection  coefficient  given  by  (2,1). 

The  infinite  iteration  of  the  system  (2,2)  starting  with  u  =  u(A),  d^^  0 

u(z)  .u(A)  ..(^)£^j^dy2^^(y^^)e-l(y^)  \  <^^,^<^i72,,^)e(y^^^^)j  .„    j  <iy^cKy^)eiy^) 
(2.U) 


'^^Sn  ^y2n-iy2 


d(z)  =  -(A)gjdy^^.,cr(y2^^^)e(y2^^)|  dy^^_2cKy2^_^)e-l(y2^2)j...  j  dy,a(y^)e(y^), 
^  ^2n-l  y2n-2  ^2 


-  u  - 

/A 
where  oij)   »  ^'(y)/4^(y)  and  e(y)  =  exp[2j  k(s)ds].  To  compute  the  reflection 

coefficient  from  (2,U)  we  must  first  set  z  =  0  and  take  into  account  the  fact 
that  o<y)  -  0  for  all  y  >  5  since  Ji(y)  is  constant  when  y  >  6.  Thus,  the  outside 
integrals  go  from  0  to  5,  and  the  inside  integrals  all  have  5  for  their  upper 
limits. 

Using  the  fact  that 

y 

(2.5)  J  or(s)ds  =  log[^(y)/^(0)]   , 

o 

we  can  integrate  the  integrals  in  (2.U)  by  parts  starting  from  the  outer 
integral  in  each  term.  Thus  we  have,  using  (2.5), 

J  -  j  dy^or(y^)e(y^)(  dy„.ia<y^_^)e-^(y^_^)(. ..  j  dy^a(y^)e(y^) 

-  e(y^)log[^(y^)/{(0)]  .  j  dy^iO^(y^,^)e"^(y^_^)j...   dy^a(y^)e(y^)l  ^  ^ 
6  6    6 

r 


+ 

L 


^yr>.i^^yn-i^^°s&^yn-i^^(°^]  *  ^^yn-i^^'^^yn-i^j—  j  ^^i^^yi^^^yi^ 

^n-l  ^2 


+  0(6), 

where  the  first  term  on  the  right  of  the  equation  for  J  vanishes.  The  term 
0(6)  refers  to  the  integral  of  a  bounded  function  between  the  limits  0  and  6. 
The  integral  J  can  then  be  written: 
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6     6 


J  -  jdyn.lO^(yn-l^^°sfi(yn.l)^(0)]   J  ...  jdy3^a(y^)e(y^)  +  0(6) 
o  y_  ,  y« 


^n-1  ^2 


6  /    / 

-  (l/2)Jdy^.ld/log^[^(y^^^)^(0)]   )/4y|  ...  (dy3_cr(y^)e(y^)+0(6). 

We  integrate  the  outer  integral  by  parts,  again  applying  the  same  sort  of 
discussion,  and  obtain 

6  6    5 

k-2 
Continuing  this  process,  we  finally  arrive  at  the  result 

1/2 

(2.6)  J  «  (l/nJ)log"[i;(6)/^(0)]   +  0(6). 

The  term  0(6)  in  (2,6)  is  the  sum  of  integrals  of  bounded  functions  between 
0  and  6,  and  each  term  of  this  sum  is  also  of  order  l/(n-l)l  by  arguments 
similar  to  those  just  used  here  and  in  report  WP-1. 

We  now  substitute  (2,6)  into  (2,1^),  remembering  our  definition  of  J, 
and  obtain 

(2.7)  u(6)  «  u(A)  +  uCA)  r   [l/(2Q)l]log^°[^(6)/^(0)]    +  0(6) 


d(6)  =  u(ZvJ  r  [l/(2n-l)!]log2^^[^(5)/^(0)]         +  0(6), 
where  we  have  used  the  fact  that  ^(y)  is  constant  for  y  >  5, 
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Tlie  quantity  ^(O)  of  equations  (2.7)  is  equal  to  ^^   of  equation  (2.1) 
and  the  quantity  5^(5)  of  (2.7)  is  equal  to  X.^  o^  (2.1).     The  term  0(5)  occurring 
in  both  of  equations  (2,7)  approaches  zero  in  the  limit  as  5  approaches  zero 
since  0(5)  is  the  stun  of  an  infinite  series, the  nth  term  of  which  is 
0[5/(n-l)i].     Therefore,   (2.7)  has  the  limiting  form 


(2.8)  u(0*)  -  r  [l/(2n)l]log2"[^2^J         -  2  cosh[log(^2'^^^^'^^] 


d(0*)  -  r  [l/(2n-l)j]log2"-^L^J       =  2  sinh[log(i:2^i)^/^]. 
n^l 

After  substituting  (2.8)  into  (2.3),  we  obtain 

(2.9)  R  -  sinh(^2^j^^^^^^sh^^2^1^^^^ 

-  [(^2^,)=^/^-  (^i/:2^^i/[(^2^i)'^'*  (^i/^^^'i 

-[V^l]/[V^l]' 

which  is  the  same  as  (2.1). 

The  small  variation  method  of  report  WP-1  gives  the  first  term  of  the 
cosh  series  in  (2.8)  for  u(0)  and  the  first  term  of  the  sinh  series  in  (2.8)  for 
d(0)  so  that 

(2.10)  R  --  log(^2^1^^^^- 

The  method  of  this  example  can  be  used  to  apply  the  small  variation 
method  to  the  reflection  problem  for  an  arbitrary  piecewise  continuous  medixm. 
Such  a  problem  must  again  be  considered  to  be  the  limiting  case  of  a  sequence 
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of  problems  for  which  the  medivun  is  continuous.     Thus,  suppose  as  usual  that 
^(z)  is  constant  for  z  >A,  and  suppose  that  J^Cz)  has  a  finite  jump  discontinuity 
at  each  of  the  points  z. , 

■i  "^  ^i+l»      0  <  z^  <A  ,  i  =  1,  ...,  n, 

where  n  is  the  number  of  discontinuities  of  5i(z)  in  the  interval  0  <  z  <  A. 
We  have  for  the  first  order  approximation  to  the  reflection  coefficient  in  the 
case  of  a  continuous  J^(z): 


(2.11) 


R--j(^'(y)/2^(y)/  exp[2jj  k(s)ds|dy 


.-6 


L  O 


n-1 


^*i*l-^  /^ 


„.+  6    z  +6  J 
i       n 


S^'(y)/4:(y)jexp(2j  j  k(s)dsjdy 

7 


-A 


^'(y)/2^(y)ye3tpl2j   k(s)dsMy. 


Now  we  consider  the  problem  P  with  the  discontinuous  5^(z)  to  be  the  limiting 
Case  of  a  sequence  of  pixjblems  Pg.  In  every  P^,  ^(z)  is  continuous  and  in 
particular  goes  continuously  from  its  value  at  z.-6  to  its  value  at  z.+6  in 
the  neighborhood  of  each  point  z.  which  is  a  point  of  discontinuity  of  the 
discontinuous  S^(z)  of  our  original  problem  P  ,  We  assume  that  the  characteristic 
impedance  ^(z)  of  each  problem  P,  of  our  sequence  is  equal  to  that  of  problem 
P  at  every  point  outside  the  intervals 


z.  -6<z<z.  +6,  i"l,  •.*,  n 
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about  the  points  of  discontinuity  z.,   i  =>  1,   ..*,   n.     We  expect  the  sequence 
of  reflection  coefficients  Re   cori^sponding  to   the  problems  P,  to  approach  the 
reflection  coefficient  R     of  problem  P     in  the  limit  as  6  approaches  zero.     To 
evaluate  this  limit  we  first  integrate  (2,11),  by  parts,  as  in  our  previous 
example .     We  have  in  the  limit  6-^0: 


1/2 


(2.12)      R/^|log  K  (Z3^))exp|2jj  k(s)dsj  -  |log  ^  (0)jexp|2jj  k(s)dsj 

=1 


n-1 


[log^  (z^^^)jexp!2j     k(s)dsj-Ws:  (z*)jexp|2jj  k(s)d£ 

^i+1  ^1 


_   ,  1/2^  1/2    )        / 

*T  (l°g  ^  (zi)  -  log  ^  (z^)}exp{2j 


i=j 


•A 
k(s)d£ 


+  jl  k(y)flog^(y)?expf2j    k(s)ds)dy. 


/ 


By  combining  terms  in  (2.12)  we  arrive  at 

(2.13)      R—  [log  tL  (z*)jexp|2jj  k(3)dsj  -  jlog  K  (0")>exp|2:)|  k(s)ds; 

z 

n 


jj  k(y)|log  i:(y)je^<2j(  k(s)dsydy. 


The  formula  (2.13)  is  completely  general  for  media  with  piecewise  con- 
tinuous Si(z).  Since  no  derivatives  appear  in  (2.13)>  this  result  can  be  applied 
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directly  to  physical  data  without  making  ary  assumptions  about  the  derivative  of 
^(z)  which  cannot  be  measured  in  practice. 

3«  Perturbation  about  a  standard  transmission  line 

Since  there  are  occasions  when  the  error  in  the  small  variation  method 
of  report  WP-1  is  not  too  large,  it  may  be  useful  to  consider  another  more 
general  perturbation  method.     If  we  can  calculate  the  wave  solutions  of  some 
standard  transmission  line  whose  electrical  properties  do  not  differ  much  from 
the  one  in  which  we  are  actually  interested,  we  can  devise  a  perturbation  scheme 
for  obtaining  approximate  wave  solutions  for  the  actual  transmission  line  of 
interest  in  terms  of  the  solutions  for  the  standard  line.     A  convenient  choice 
for  a  standaixi  transmission  line  from  which  we  can  construct  a  perturbation 
solution  is,   for  example,   a  homogeneous  line  whose  constant  electirical  parajtneters 
L  and  C  are  perhaps  average  values  of  the  L(z)  and  C(z)  of  the  actual  line 
which  may  be  inhomogeneous.     Some  of  the  advantages  of  the  small  variation 
method  will  not  be  present  in  the  new  method;  e.g.,  the  error  estimate  of  our 
former  approximation  was  independent  of  frequency,   as  we  observed  in  report 
WP-1,  whereas  the  error  in  the  present  approximation  will  depend  upon  the 
frequency  in  general.     The  present  method  can  then  be  expected  to  provide  a 
more  appropriate  solution  to  the  inhomogeneous  transmission  line  problem,  only 
in  the  highly  specialized  situation  for  which  it  is  designed  in  a  given  case. 

We  start,   as  before,  with  the  problem  of  obtaining  the  wave  solutions 
V(z)  and  I(z)  of  an  inhomogeneous  transmission  line  system 

(3.1)  d7(z)^z  =  -  j«I(z) 

dl(z)Az  =  -  ja)C(z)V(z), 

where  we  have  chosen  units  of  length  so  that  L(z)  =■  1  in  (3.1).     At  some  point 
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z  "A,   a  load  impedance  Z  or  load  admittance  T  =  l/Z  is  given  so  that 

(3.2)  1(A)  =  TV  (A). 

Let  US  assume  that  we  have  available  the  complete  set  of  linearly 
independent  solutions  of  a  system 

(3.3)  dv(z)/dz  =  -  >i(z) 

di(z)/dz  =  -  Ja)/^(z)i(z)} 

i.e.,  we  have  a  non-singular  solution  matrix 

(3.U)  /v^(z)        V2(z)\ 

which  satisfies 

(3.5)  d(v)/dz  -  -  j«(K)(v), 

(o     l] 

where  (K)  =W(z)  O/. 

It  is  a  well-known  fact  that  the  matrix  (v)  of  (3.'i)  has  a  constant 
determinant  for  all  z.     This  determinant  is  fixed  arbitrarily  by  properly 
choosing  the  initial  values  of  the  independent  solutions  which  form  its  columns, 
We  can  therefore  require  that 

(3.6)  det(v)  =  1. 

Relation  (3.6)  implies  that  the  inverse  matrix  of  (3. '4)  is  given  by 

(3.7)  T     /i,   -V, 
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If  we  let  E  be  the  vector  [  ,  j,  (C)  the  matrix  (_  q|  and  (P)   the 
matrix  (C)  -  (K),  we  can  write  (3«1)  in  the  form 

(3.8)  dE/dz  -  -  3co(C)E  =  -  XK)E  -  XDE. 

The  equation  (3.8)  has  a  form  suitable  for  applying  the  perturbation 
method,  perturbing  about  the  standard  transmission  line  of  (3.5).  To  carry 
out  the  process  we  write 

(3.9)  E  -  (v)E, 
where,  for  the  components  of  E,  we  have 


On  substituting  (3.9)  into  (3.8),  making  use  of  (3«5)»  we  obtain 
(3.10)  dE^z  -  -  j<o(v)"^(r)(v)E. 

In  order  to  complete  the  perturbation  procedure  we  shall  form  an 

A 

equivalent  system  of  integral  equations  for  the  components  of  E.  However,  to 
obtain  the  most  rapidly  convergent  iteration  sequence  we  can,  without  complicating 
the  final  expressions,  we  should  remove  the  diagonal  part  of  the  matrix 
(▼)~  (T')(v)  in  (3.10)  and  further  normalize  the  still  somewhat  arbitrary  matrix 
(v)  so  that  one  component  of  the  initial  vector  in  our  final  system  of  integral 
equations  is  zero.  As  we  have  observed  in  Section  2,  and  in  report  WP-1,  a 
system  of  integral  equations,  with  an  initial  vector  having  one  component  zero 
and  an  anti-diagonal  kernel  matrix,  will  produce  iteration  series  for  the 
components  of  the  solution  vector  such  that  each  series  will  be  lacking  every 
other  term. 
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With  this  principle  in  mind,  we  take  (a)  the  diagonal  part  of  the 
matrix  (v)"  CP)(v)  of  (3.10),  and  define  a  diagonal  matrix  (7-)  as  a  solution 
of  the  differential  equation 

(3.11)  d(r)/dz  =  -  M^)M' 
We  then  define  E  by 

(3.12)  E  =  (r)E. 

After  substituting  (3.12)  into  (3.10)  we  obtain,  using  (3.11) 

(3.13)  dE/dz  =  -  Xr)'-'-(Y)(r)E, 

where  (y)  is  (v)~  (7^)(v)-(cr),  and  we  can  write  finally: 

dE/dz  =  -  jco(G)E, 

where   (G)  =  ('?-)"  (y)(T')  is  an  anti-diagonal  matrix,   as  required. 

If  we  choose  initial  conditions  on  the  matrix  (v)  of   (3.U)  in  such  a 
way  that,   in  addition  to  det  (v)   =  1,  we  have 

(3.11;)  i-^CA)  =  TVj^Ca), 

the  component  I(Z^  of  the  vector  E  at  z  =  A  will  vanish;  i.e.,  we  shall  have 

(3.15)  1(A)  .(^^'\. 

as  required. 

Before  carrying  out  the  last  step  in  deriving  our  approximate  wave 
solution  for  the  transmission  line  given  by  (3.1) *  let  us  set  down  explicitly 
the  expressions  for  the  various  matrices  which  occur  in  this  derivation. 
We  have,  writing 
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(3.16)  C(z)  -X(z)  =  Q(z), 
the  matrices: 

(3.17)  (f)-     (l^^^      I J    , 


(3.18)  (cr)  =  Q(z)v^(z)v2(z)    ("J      J 


/«q)y-jco    Q(y)vT(y)v-(y)dy  0 

(3.19)  ir)  '{       ^      i  (   f 

0  exp|>    Q(y)v^(y)v2(y)dy| 

z 


0  -v^(z)exp|2jco  v-(y)r,(y)Q(y)dyj 

(3.20)  (G)  =  Q(z)  I  ,  /2^     ^  ^      i,  ^ 


r^(z)e:5'<'-  2jto  v^(y)v2(y)Q(y)dyl  0 


After  an  integration  (3.13)  becomes 

(3.21)  E(z)  =  E(A)  +  ja>j(G)E(y)dy. 

This  leads  to   the  first  order  approximation  obtained  by  a  single  iteration, 

(3.22)  V(z)  ^  ?(A) 

l(z)^l{A)ia>\  Q(y)v^(y)exp|-2jco     Q(s)v^(s)v2(s)ds|dy. 
z  y 


.  lU  - 

The  error  of  the  approximation  (3.22)  can  be  estimated  by  the  expressions 
(U.12)  and  (U.U)  of  report  WP-1,  The  function  ^.(z),  given  by  (U.U),  is  in  the 
present  case 

(3.23)        tx(2)  =  «|Q(z)|[k^(z)|^  *   |v2(z)|HV^ 

providing  the  phase   I  "<^(y)v_  (y)Q(y)dy  is  real.     The  error  given  by  (U.12)  of 
report  WP-1  for  V(z)  is 

(3.2li)  e^  <  exp  jti(y)dy  -  Jl+L(y)dy| 

z  z 

and  for  l(z) 

^A  12 


ej  <  exp  In  (y)dy  -  f  1  +  |  n  (y)dy  + 


M.  (y)dy 
•-'z 


'2   . 


These  expressions  follow  from  the  fact  that  each  of  the  two  iteration  series 
skips  every  other  term  so  that  the  next  term  of  the  series  for  V(z)  is  of  second 
order  aiid  the  next  term  of  the  series  for  I(z)   is  of  third  order. 

U»     Example: 

As  an  example  of  the  perturbation  method  described  in  Section  3  we 
shall  consider  the  case  which  arises  when  a  plane  wave  is  incident  at  nearly 
grazing  incidence  on  a  stratified  medium.     It  is  convenient  in  this  case  to 
perturb  about  a  medium  with  C(z)  =  0,     Thus , /if  (z)  -  0,  and  we  have  for  the 
equations  (3*2) 

(U.l)  dv/dz  =  -;5a)i 

di/dz  =  0 . 
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The  general  solution  of  (U.l)  has  the  form 

(U.2)  i  =  constant 

V  =  -jcoiz  +  y(0). 

We  niust  choose  appropriate  values  of  i  and  v(0)  to  form  a  solution 
matrix  whose  determinant  has   the  value  one  and  so  that  (3.1ii)  holds.     Moreover, 
in  order  to  minimize  the  error,   the  remaining  constants  should  be  chosen  so 
that  the  phase  of  the  integrand  of  (3.22)   is  realj   i.e.,  the  product  r^(z)vAz) 
should  be  real,     A  choice  of  constants  which  satisfies  these  requirements 
gives  for  the  elements  of  the  solution  matrix 

T^(z)         VgCz) 

of  (3,3)  the  values 

(U.3)  v^(z)  =  xll  -  «Va   (z-A)  -  jcoTzV  /(I  +coVa^) 

VgCz)  -  -  jl-  A^A  (z-A)  +  jcoYzV    /2XT 

i^(z)  «  XT  /(1+ jcoTA)  (  =  constant) 

igCz)  =  (l  +  jcoT/i)   /2X  (  »  constant), 

where  the  constant  X  is  still  arbitrary. 

Since  X  is  arbitrary  it  can  be  chosen  so  that  the  quantity  n(z)  of 
(3»23)  which  appears  in  the  error  estimate  (3.2U)  is  a  minimum.     We  observe 
that  except  for  constant  factors  m.(z)  has  the  form 

(U.U)  h(z)  =  (aX^+px"^)V^ 
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from  the  form  of  the  expressions  for  -rAz)   and  t^{z)   in  (U.3).  By  differentiating 
(U.U)  with  respect  to  X  and  setting  the  resulting  expression  equal  to  zero,  we 
obtain  for  the  minimizing  value  of  X, 

We  compare  (U.li)  with  (3.23)>  taking  into  account  (U.3),  and  obtain 
from  (1.5) 

(U.6)  X  =  (1  *  coVA2)^/y  (2Y)^/^ 

Observe  that  fortuitously  the  minlMzing  X  of  (U.6)  is  independent  of  z,  as 
it  must  be  if  it  is  to  be  of  any  use,  since  we  can  only  allow  a  constant  X 

in  (U.3). 

In  the  present  example  the  function  Q(z)  is  given  by 

(U.7)  Q(z)  -  C(z). 

From  the  following  expressions  the  approximate  wave  solutions  of  (3.1), 
subject  to  the  load  admittance  T  at  the  point  z  =A,  can  be  obtained: 

(U.8)        V(z)~'  V(A) 

i(z)^  V(a)  |C(y)v^(y)exp  2jcojc(s)|v^(s)lWdy, 

where 

(U.9)  V3_(z)  -  |l-  co2y2a(z  -a)  -  jcoTzj /|2T(1  +a)VA^)J-^^. 

To  obtain  V(z)  and  l(z)  we  combine  (3.9)  and  (3.12): 

E  -  (v)(t-)E. 
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We  then  have 
(li.lO) 

V(z)'^V(A) 


/A 


r^  (  z)eiq?  j"  jco    C(y)  k^^Cy)  |  dy 


-3'^v^(z)Jc(y)v^(y)expi2jco    C(s)|v^(s)|^ds|dy  exp^ 


z 


-  jcojc(y)|v^(y)|V} 


I(z)~V(A)(T/2) 


1/2 


e  expJjco|c(y)|v^(y)|dy 


+  aoje|c(y)v^(y)  exp/2jco  jc  (3)  !v^(s)  l^dsidy  asp/- jw  j  C(y)  Iv^(y)  |dy> 

where  9  =  (1  +  coTf  A  )        /  (1+ j«YA),  a  complex  number  of  unit  modxilus. 

If  it  is  desired,  the  reflection  coefficient  can  be  computed  from  the 
relation 


(U.n) 


R  =  [v(o)c^/^(o)  -  i(o)]/[v(o)c^/2(o)  +  1(0)]. 


As  an  example  of  how  one  might  use  the  formiilas  to  estimate  quantities 

of  interest,  consider  the  case  of  a  very  thin  layer:  A '^  0.  From  (U.ID), 

2 
(U.3)  and  (U,6)  we  get  roughly  (neglecting  terms  of  order  A  ): 

(U.12)        V(0)^  V(a)  /(2Y)^/^ 

1(0)^  V(a)[y  +  jcoa(c(A)  -t2]]/(2Y)^/2. 

2 
If  we  assume,  as  is  the  usual  case,  that  C(A)  =  T  ,  then  we  have 

(U0I3)         R  --  [c^/^  (0)  -  T  ]  /  [c^/^  (0)  +  y]  +  O(A^). 
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Except  for  terms  of  order  A  >   (I4..I3)  is  just  the  expression  for  the  reflection 
coefficient  of  two  different  homogeneous  trsnsmission  lines  joined  at  the 
point  z  =  0,  where  1/c  '    (O)  is  the  characteristic  impedance  of  the  fii^t  line 
and  lA  is  the  characteristic  impedance  of  the  second  line.     In  the  particular 
case  where 

C(0)   =  C(A), 

2 

R  vanishes  except  for  terms  of  order  A  . 
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